Expansion of a quantum gas released from an optical lattice 
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We analyze the interference pattern produced by ultracold atoms released from an optical lattice. 
Such interference patterns are commonly interpreted as the momentum distributions of the trapped 
quantum gas. We show that for finite time-of-flights the resulting density distribution can, however, 
be significantly altered, similar to a near-field diffraction regime in optics. We illustrate our findings 
with a simple model and realistic quantum Monte Carlo simulations for bosonic atoms, and compare 
the latter to experiments. 
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Experiments with ultracold quantum gases in optical 
lattices rely heavily on time-of-flight (ToF) expansion to 
probe the spatial coherence properties of the trapped gas 
E B B 0, 0, S, 0, H ■ When the phase coherence length is 
large compared to the lattice spacing, the post-expansion 
density distribution shows a sharp interference pattern 
with the same symmetry as the reciprocal lattice. As 
the phase coherence length decreases, e.g., on approach- 
ing the Mott insulator (MI) transition, the visibility of 
this interference pattern decreases accordingly [lj]. To 
obtain a more precise understanding beyond this quali- 
tative description, it is usually assumed that the density 
distribution 71tof(f) of freely expanding clouds provides 
a faithful map of the initial momentum distribution. 

In this Letter, we point out that, in general, the ToF 
distribution differs from the momentum distribution for 
finite time-of-flight, the latter being recovered only in 
the "far-field" limit t — > oo. Practically, the ToF and 
momentum distributions become identical after a char- 
acteristic expansion time tpp = mR$l c /h, which depends 
on the particle mass m, the coherence length l c , and the 
cloud size Rq prior to expansion. This time scale can be 
understood in analogy with the diffraction of a coherent 
optical wave by a periodic grating. Then, the character- 
istic £ff in the expansion problem exactly corresponds 
to the Fresnel distance in the diffraction problem. The 
far-field regime is typically reached when the coherence 
length is short, for example for a cloud in the MI regime, 
or a thermal gas well above the critical temperature. We 
show that for phase-coherent samples where a sizeable 
fraction of the atoms are Bose condensed, the far-field 
condition is usually not met for typical expansion times 
used in current experiments [H,[lS,S,S, 0,0,1- Exper- 
imental measurements and quantum Monte-Carlo simu- 
lations are used to demonstrate that this results in sub- 



stantial changes in the ToF distribution. We also discuss 
implications for the interpretation of the ToF images. 

We consider an ultracold boson cloud released from 
a periodic trapping potential with cubic symmetry, lat- 
tice spacing d = Al/2, and lattice depth Vq given in 
units of the single-photon recoil energy £r = /i 2 /2mA 2 , 
where Al is the lattice laser wavelength. In addition to 
the lattice potential, an "external" harmonic potential 
is present, due to both the magnetic trap and the opti- 
cal confinement provided by the Gaussian-shaped lattice 
beams 0, H[ • This external potential is responsible for 
the appearance of a shell structure of alternating MI and 
superfluid regions in the strongly interacting regime. 

The density distribution after expansion for a time t is 
usually expressed as a product (see, e.g., [Toj] ) . 
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where an envelope function wq is the Fourier transform 
of the on-site Wannier function wq) and the interference 
term is 



5(k) = £ e 
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Here the operator a J, creates an atom at site r M . To assess 
the validity of the far- field approximation used in Eq. |T]), 
we quickly outline its derivation. Neglecting interactions 
during expansion (see below), the atomic field operator 
can be expressed in Schrodinger's picture as ^(r,t) = 
Y,r v W v (r,t)a v where W„(r,t = 0) = w (r - r„). After 
the cloud is released, the wavefunction W v evolves in free 
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LOht^> 1, with ujl the oscillation frequency at the bottom 
of a lattice well. In the limit t — > oo, the dependence 
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FIG. 1: (a) Momentum distributions for a one-dimensional 
lattice with parabolic distribution of the occupation numbers 
calculated using Eq. ([5]) (solid line, expansion time t = 20 ms, 
dashed line: expansion time t — 100 ms, dot-dashed line: 
expansion time t — > oo). (b) Evolution of the peak amplitude 
A with expansion time t/tpF. The dashed line shows the 
expected near-field scaling in one dimension, A oc t/t-pF- The 
number of sites is 2JVtf+1 = 61 for (a) and (b). (c) Evolution 
of the width of the diffraction peaks with expansion time. 
The width has been normalized to the separation between 
two adjacent diffraction peaks for convenience. The circles 
show the experimental measurements and the solid line a fit 
by a hyperbola oc 1/t, as expected in the near-field. 



on the initial site position r„ vanishes, and one recovers 
Eq. |T|). For finite t, this dependence can be neglected in 
the envelope function ll|, but not in the phase factor. 
We thus obtain a generalized interference term 



(3) 



Note that experimentally one observes a column distri- 
bution integrated along the probe direction, Sj_(k^) = 
/ dk z \w(k z )\ 2 St(k). This is included in latter compari- 
son with experiments, but in the following we base our 
discussion on Eq. ^ for simplicity. 

A fruitful analogy can be made with the theory of opti- 
cal diffraction. The formation of the interference pattern 
results from the interference of many spherical matter 
waves emitted from each lattice site, with phase relation- 
ships reflecting the initial quantum state of the boson 
gas. We can exploit this analogy further by defining the 
equivalent of a Fresnel distance usually introduced in the 
theory of optical diffraction to estimate the importance of 
the quadratic phase factor oc — r 2 . Because the corre- 
lation function (fl* a v ) suppresses contributions from sites 
distant by more than the characteristic coherence length 
l c , we can estimate the magnitude of the quadratic phase 



and ~ m fi/*° near the cloud edge. Here Rq the char- 
acteristic size of the cloud before expansion. The most 
restrictive condition to apply the far-field approximation 
thus reads t t FF , with 
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As an example, for a 87 Rb condensate with l c m Rq ~ 
30 d and a lattice spacing d ~ 400 nm, one finds t FF ~ 
100 ms, much larger than typical expansion times t ~ 20 
ms in experiments 12j. In contrast, a gas with short co- 
herence length (e.g., in the MI regime), with l c > d, will 
enter the far-field regime after a few ms. We stress that 
the quadratic Fresnel term is intrinsically non-local, as 
the dephasing between two particular points r M and r„ 
depends not only on their relative separation but also 
on their absolute positions. Although this has little ef- 
fect deep in the superfluid or in the MI phase, this casts 
serious doubts on the validity of a local density approx- 
imation to compute quantitatively the ToF distribution 
in regimes where the coherence length is intermediate be- 
tween the cloud radius and the lattice spacing. 

To illustrate the influence of Fresnel terms on the in- 
terference pattern, we consider a ID lattice with uniform 
phase and parabolic distribution of the occupation num- 
bers, (fijjfij/) = c^Cy with c M = v/l — (/i/iVxp) 2 - The ToF 
distribution is given by 
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with k = kd, Ntf = Ro/d — 30 the Thomas- Fermi con- 
densate size in lattice units, and where (3 = y/md 2 /ht. 
The normalization factor (2Ntf + l) 2 would give the 
peak amplitude if the filling factor were uniform. We 
plot in Fig. the distributions corresponding to t = 
20 ms~ 0.2iFF which shows a significant broadening of 
the distribution for short time of flight when compared 
to the asymptotic result. For longer expansion times 
t ~ £ff ~ 100 ms, the far-field approximation is recov- 
ered to a good approximation. 

Qualitatively, we expect from dimensional arguments 
that the peak width scales as (f3N TF ) 2 = t FF /t in the 
near-field, while approaching a constant value in the far- 
field. The peak height thus increases as {t/t-p F ) T> in T> 
dimensions. This is confirmed by the one-dimensional 
calculation shown in Fig. [TJd. This dependence provides 
a mean to check the importance of near-field effects ex- 
perimentally. For the measurement, a sample of roughly 
10 5 87 Rb atoms has been been prepared in a three- 
dimensional optical lattice with a de pth V y — 6 Er, and 
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subsequently released for expansion 14| . After record- 
ing a series of absorption images for different expansion 
times, the width of the interference peaks was extracted 
using a Gaussian fit to the images. We plot the results 
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in Fig. [TJ;, normalized to the separation between two 
diffraction peaks for convenience. The data confirms the 
t FF /t scaling, indicating that the far-field asymptote is 
not reached even after the longest expansion time avail- 
able in the experiment. 

We now discuss briefly the effect of interactions on the 
expansion, and show that this is negligible compared to 
the finite ToF effect. When the cloud has just been re- 
leased from the lattice potential, each on-site wavefunc- 
tion W/j, expands independently with a characteristic ex- 
pansion time lo^ 1 , until t f» t* = yjh/ (lu^Er) where 
the wavefunctions expanding from neighboring sites start 
to overlap. At this time, in the usual situation where 
ujjJj* ^> 1, the local density has dropped dramatically by 
a factor (wj,i) -3 <C 1. Hence, the interaction energy con- 
verts into kinetic energy on the time scale of a few oscilla- 
tion periods only, and expansion becomes rapidly ballis- 
tic. The parameter controlling the importance of interac- 

tions is given by V = fa « y/te eg* (H)^ with U 
being the on-site interaction energy. For typical parame- 
ters, T) is small (for instance rj ~ 0.05 for Vq = 10 En and 
the experimental parameters of Q). Hence, we expect 
only small corrections to the non-interacting picture of 
ballistic expansion. This has been confirmed using a vari- 
ational model of the expanding condensate wavefunction 
[l~5| . This model predicts that the "Wannier" envelope 
expands faster as compared to the non-interacting case, 
which does not affect the interference pattern, and picks 
up a site-dependent phase factor formally similar to the 
Fresnel term discussed previously, but with a very weak 
prefactor rj <C 1 which has negligible influence in prac- 
tice. We conclude that interactions essentially contribute 
to the expansion of the on-site wavefunctions, without 
significant dephasing of the interference pattern. 

The discussion so far focused on fully phase-coherent 
systems, which only applies to the weakly-interacting 
regime at low lattice depths. To investigate how the 
interference pattern is affected for strongly interacting 
systems (i.e., on approaching the Mott transition and be- 
yond), we have performed large-scale three-dimensional 
quantum Monte Carlo (QMC) simulations accounting 
for the external trapping potential using the worm al- 
gorithm 1(3, 12 1 m the implementation of Ref. fl8| . The 



calculations were performed for N = 8 x 10 4 atoms, us- 
ing exactly the same parameters and system sizes (up to 
~ 200 3 ) as in the experiments reported in The simu- 
lation was done at low constant temperature T = J/ks, 
where J is the hopping amplitude. Although simulations 
at constant entropy would be closer to the experimental 
situation, the temperature turns out to be approximately 
constant in this parameter regime [l9j |. 

The ToF distribution calculated for finite and infinite 
expansion times are shown in Fig. [2] The simulations 
confirm explicitly the analysis made above: the interfer- 
ence pattern is strongly affected in the superfluid phase, 
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FIG. 2: Results from Quantum Monte Carlo simulations. On 
the left column, we show a horizontal cut through the ToF 
distributions for a finite expansion time t — 14 ms (solid line), 
compared to a cut through the profile calculated for t — > do 
(dashed line). Units for n± are arbitrary. The insets show 
directly the two-dimensional ToF distributions for t = 14 ms. 
On the right column, we show the in-trap density profiles for 
reference. The lattice depths are Vb = 12-Er (&)d), 15Er 
(b,e) and 17Er (c,f), respectively. 



and the effect becomes less and less pronounced as the 
lattice depth is increased and the Mott transition crossed. 
Note finally that the Fresnel phase suppresses the contri- 
bution from the edges of the cloud, thus favoring the con- 
tribution of the central region to the ToF pattern. This 
is especially important when superfluid rings surround a 
central MI region with lower coherence [201 ] . 

The interference p_attern is often characterized by its 
visibility 0,11 US 



V 



"ToF(k max ) ~ ™ToF(k m i n ) 
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(2-7T, 0) and k m i n e? = V2(n, it) to 
cancel out the Wannier envelope in the division. We first 
evaluate the sensitivity of V to the Fresnel phase by plot- 
ting in Fig. ((3]) two theoretical "benchmark" curves as- 
suming perfect experimental resolution (dashed and dot- 
dashed lines for t = 14 ms and t — ► oo, respectively). We 
find little difference between the two curves when T j J is 
kept constant and small. Indeed, the Fresnel terms only 
matter for systems with large coherence length, where the 
visibility is by construction very close to unity. We con- 
clude that a detailed investigation of the superfluid side 
of the transition is better achieved by directly measuring 
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Lattice depth V (f ft ) 



This could eventually provide a way to investigate the 
physics near the quantum-critical point without "par- 
asitic" contributions coming from coexisting superfluid 
rings. 

Simulations were ran on the Brutus cluster at ETH 
Zurich. We acknowledge support from IFRAF, ANR 
(FG), DFG, EU, AFOSR (IB), the Swiss National 
Science Foundation (LP), NSF grant PHY-0653183 
(BCS,NP,BS) and DARPA (OLE project). 



FIG. 3: Visibility of the interference pattern as defined in 
Eq. (|6]). The dashed and dot-dashed lines show the Quan- 
tum Monte Carlo result for infinite and finite (t — 14 ms) 
expansion times, assuming perfect experimental resolution. 
The solid line is computed for t — 14 ms accounting for finite 
experimental resolution. Note that the comparison between 
experiment and simulation is only qualitative, since the sim- 
ulations were performed at constant temperature T = J/Ub 
while the experiment was not. 



the ToF distributions, whereas the visibility is well-suited 
for short coherence lengths. 

We also compare in Fig. [3] the experiments reported in 
Q to the predictions of the QMC simulations (solid line). 
Here, we emphasize that apart from the Fresnel terms, an 
accurate comparison requires to account for the experi- 
mental resolution, which is limited by two effects. First, 
the signal was obtained by integration over a square box 
centered around the maxima or minima, the integration 
area being ps (0.11 x 2ir/d) 2 in momentum units. This is 
comparable to a typical peak area, so that the visibility 
is calculated from the peak weight, rather than from its 
amplitude. Second, the finite resolution of the imaging 
system (about 6 /im) is not negligible for the sharpest 
peaks. Accounting for these two effects when evaluating 
the QMC data, we find good agreement with the exper- 
imental results. This entails that the experimental data 
arc compatible with the system remaining at low enough 
temperatures to cross a quantum-critical regime, in con- 
trast to the analysis made in Rcfs. 



21. 22] which included 



neither near-field expansion nor experimental resolution. 

In conclusion, we have analyzed the interference pat- 
tern observed in the expansion of a bosonic quantum gas 
released from an optical lattice. We showed that due to 
an additional Fresnel-like phase appearing for finite time 
of flight, the ToF distribution can be markedly different 
from the momentum distribution for clouds with large co- 
herence lengths. Conversely, the visibility as calculated 
from Eq. ^ is rather insensitive to this effect. 

The Fresnel phase acts as a magnifying lens for the 
central region undergoing a Mott insulator transition by 
suppressing the contribution of the outer regions of the 
cloud when the central density is close to integer filling. 
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